Recombination limited energy relaxation in a BCS superconductor 
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We study quasiparticle energy relaxation at sub-kelvin temperatures by injecting hot electrons 
into an aluminium island and measuring the energy flux from electrons into phonons both in the 
superconducting and in the normal state. The data show strong reduction of the flux at low tem- 
peratures in the superconducting state, in qualitative agreement with the presented quasiclassical 
theory for clean superconductors. Quantitatively, the energy flux exceeds that from the theory 
both in the superconducting and in the normal state, possibly suggesting an enhanced or additional 
relaxation process. 
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Superconducting nanostructures attract lots of atten- 
tion currently, partly because of their potential appli- 
cations, for instance, in single Cooper pair and single- 
electron devices, in quantum information processing, and 
in detection of radiation. Although the operation of 
many of these devices is based on charge transport, the 
energy relaxation in them is also of importance to war- 
rant proper functioning either under driven conditions, or 
when subjected to environment fluctuations. Thermal- 
ization of the electron system with the surrounding bath 
is a serious concern at sub-kelvin temperatures for non- 
superconducting structures, but securing proper thermal- 
ization of a superconductor is an even greater challenge. 
In particular, recombination of hot quasiparticles (QP:s) 
into Cooper pairs slows down exponentially towards low 
temperatures. QP scattering rates in usual BCS super- 
conductors have been assessed theoretically already sev- 
eral decades ago [1,, t2j], and there have been measure- 
ments of them, both soon after the first predictions (for 
review, see ll| ) , and recently also at very low tempera- 
tures d, 0, U However, although a well recognized issue 
in normal systems, the most relevant property of relax- 
ation, the associated heat flux, has not been addressed 
in the past. This is the topic of the present letter. We 
present both experimental and theoretical results which 
demonstrate the importance of slow thermal relaxation 
in superconducting nanostructures. 

Energy relaxation in normal metals has been investi- 
gated thoroughly in experiment for a long time @, 0,11,0] ■ 
The central results can be summarized as follows. In 
three dimensional systems electron-phonon (c-p) heat 



flux Pep is 



Pep = SV(T5 - 
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Here, E is a material constant V is the volume of the 
electronic system, and Pc and Pp are the temperatures 
of electrons and phonons, respectively. Deviations from 
this behaviour towards the fourth power of temperature 
have been seen for lower temperatures (see, for example. 




FIG. 1: A typical sample (Sample C) for measuring energy 
relaxation in an aluminium superconducting bar. The circuits 
on the right indicate injection of hot QP:s and probing the 
island temperature. 



llll), and are usu- 



Ref. [ij and also more recent Ref. 
ally explained by the impurity effects [l^, |13|, when 
the wave length of a thermal phonon becomes of the or- 
der of the electron mean free path or of the sample size. 
Nevertheless, Eq. ([T]) gives a good account of the heat 
flux observed in most experiments at sub-kelvin temper- 
atures. Under the same conditions electron-electron (e- 
e) relaxation is typically much faster; most experiments 
demonstrate the so-called quasi-equilibrium, where elec- 
trons have a well-defined temperature, usually different 
from that of the phonons. Deviations from this simple 
picture have been observed, e.g., in voltage biased diffu- 
sive wires [l5| . 

Relaxation processes in superconductors have also 
been studied, see, e.g., Ref. The most obvious fea- 
tures different from the normal state are: (i) The QP:s 
need to emit or absorb an energy in excess of the energy 
gap A to be recombined or excited, respectively. This 
leads to exponentially slow e-p relaxation rates at low 
temperatures, (ii) The number of QP:s is very small well 
below Tc, leading to slow e-e relaxation as well. The 
focus has been in relaxation times, with no attempts to 
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obtain energy flux in the spirit of Eq. ^ . The relaxation 
time was addressed recently, e.g., in experiments on su- 
perconducting photon detectors [j, 0, [^; these measure- 
ments suggest to confirm the recombination limited rate 
oc y^T/Tce-^/'"'^ down to T/Tc ~ 0.2. At lower 
T the relaxation time saturates due to presently poorly 
known reasons. 

The e-p processes in clean superconductors can be 
characterized by the rate Tj^^-q '^^ ^ Q-P with wave 
vector k to emit a phonon with wave vector q, 
T-k.l-ci = (27r/n)|A^k.k-qP<5(i^k-£^k-q-eq)iV(i;k-q)[l- 
f(E]^-q)][np{q,Tp) + 1]. Here, >lk,k--q is the ma- 
trix element of electron-phonon coupling in a supercon- 
ductor, f{E) is the distribution function of electrons; 
it is the Fermi function f{E,T^) = (1 + e^/'^^^-)'! 
if electrons are in equilibrium with a temperature T^. 
Phonons are assumed to be in equilibrium with oc- 
cupation np(q, Tp) — (^e'^i/'^sTp _ 1)^^. Compared 
to the normal state we have inserted the normal- 
ized density of states (DOS) N{E). For a supercon- 
ductor with energy gap A(r), the DOS also depends 
on temperature N{E,T) = \E\/yjE'^ - A{T)^Q{E^ - 
A(T)^), where Q{x) is the Heaviside step func- 
tion. Electrons emit energy to phonons at the rate 
Pe = N{Ef) J dE^N{Ek)f{E^,%) J d^qDp{q)e^T^l_^. 
Dp{q) is the phonon DOS and N{Ef) is the normal-state 
DOS at the Fermi level. Writing a similar expression for 
absorption of phonons Pa, we obtain the net heat flux, 
Pep = Pe — Pa- Calculating the rate and the matrix 
elements from the quasiclassical theory for clean super- 
conductors [iBl we flnd the e-p heat flux 



TABLE I: Sample dimensions and junction resistances. 
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which is by a factor 0.98e^'^/'^^"^'' smaller than the result 
for the normal state [Eq. (P) with Tp < To]. 

Figure [T] shows a typical configuration of our experi- 
ments. The samples were made by electron beam lithog- 
raphy and shadow evaporation. The parameters of the 
structures are given in Table [H The aluminium block in 
the centre of Fig. [1] is the volume in which energy re- 
laxation is investigated. Two small and two large tunnel 
junctions connect the island to aluminium leads. The hot 
QP:s are injected via one of the small tunnel junctions 
in series with a large one. Because of the large asym- 
metry of junction parameters, essentially all the power is 
injected by the small junction. The steady-state distri- 
bution on the island is deduced from the current-voltage 
curves (IVs) of the opposite pair of junctions. We observe 



Sample 


volume (/im^) |i?i,i?2,i?3,-R4 (kf2) 


A 
B 
C 


21 ■ 1.5 • 0.44 840, 4, 4, 1160 
4.9 • 1.5 ■ 0.44 760, 5.7, 5.7, 1290 
4.9 ■ 1.5 ■ 0.44 485, 20, 20, 980 



the QP current of only the small junction; the large junc- 
tion remains in the supercurrent state. Measurements in 
a configuration with two small junctions in series as in- 
jectors and the two large junctions in series as probes 
were also made with essentially identical results. 

If the island is at temperature To and the leads at 
Toxt, the QP current / is given by eRrl = / dEN{E - 
eV,T,)N{E,T,^t)[fiE - eV,T,) - /(£;,ro,t)]. Here Rt 
is the normal state resistance of the junction. For equal 
island and lead temperatures, Te = Text, the calculated 
(a) and measured (b) IVs for various T^/Tc are shown 
in Fig. O Wide plateaus in the regime < eT^ < 2A 
emerge due to the thermal QP current; its value at 
eV = A is shown in (c). The agreement between ex- 
periment and theory is good down to Te/Tc — 0.25. 
To match the data to the theory also at lower tempera- 
tures one can use the pair-breaking parameter 7 = F/A 
which yields a smeared density of states, N{E,T) — 
\Re{E + iT)/^{E + iT)^ ~ A(To)2|. In the figure we 
show lines with 7 — 10~^ and 7 — 10~^. Since at higher 
temperatures no fit parameter is needed, we focus our 
analysis to the range 0.3 < Tc/Tc < 1 to be on the safe 
side. 

Figure[2](d) shows the calculated IVs of the probe junc- 
tion, assuming that only the island temperature Te is ele- 
vated, and the leads remain at Text — O.OSTc- This is the 
expected behaviour under power injection, provided the 
e-e relaxation is strong and that the junctions are opaque 
enough not to conduct heat from the island into the leads. 
A peak in the IVs arises at eV = A(roxt)- A(To). In Fig. 

(e) we show the corresponding measured curves at var- 
ious levels of injected power. The resemblance between 
(d) and (e) is obvious but the features in experimental 
curves are broadened in comparison to those from the 
theory, which is common for small junctions [l3|- In the 
data analysis we next find the minimum current in the 
plateau-like regime at bias voltages between the "match- 
ing" peak and the strong onset of QP current. This cur- 
rent is converted into temperature by comparing it to the 
To dependent minimum current of the theoretical IVs. 

The power Q{V) deposited on the island by a biased 
junction is given by e^RrQiV) = J{E - eV)N{E - 
eV, T,)NiE, T,^)[f{E, Text) - f[E - eV, T,)]dE. This 
equation allows us to determine the injected power, as 
well as the heat flux through all the junctions. There are 
two features to note: (i) Since typical injection voltages 
in the experiment are V S> A/e, it is sufficient to assume 
that the power injected into the island equals IV/2, i.e., 
it is divided evenly between the two sides of the junc- 
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FIG. 2: Tunnel currents under equilibrium and quasi- 
equilibrium conditions for a superconductor. Theoretical (a) 
and experimental (b) IVs of a junction at several bath temper- 
atures, (c) Theoretical and experimental currents at eV = A. 
The two theory lines correspond to two different realistic pair 
breaking parameters 7 = 10""^ (upper curve) and 7 = 10"* 
(lower curve), (d) Calculated IVs when the two leads of the 
junction have different temperatures, (e) The measured IVs 
under a few injection conditions, (f) The current in Sample 
A on the plateau between the initial peak and the rise of the 
current at the conduction threshold around 2A/e. 



tion: the junction behaves essentially as a normal junc- 
tion, where this statement is true always, (ii) The heat 
conductance of the (probing) junction is almost constant 
over a wide range of voltages within the gap region. Its 
value is low and can be neglected under most experimen- 
tal conditions. Yet, to test this, we varied the resistances 
of the large tunnel junctions by a factor of five between 
samples A and C, without a significant effect on the re- 
sults. Figure [2] (f) shows the current on the plateau, 
as described in the previous paragraph, as a function of 
power injected, at various bath temperatures. In a wide 
range, from 30 mK up to 380 mK, the behaviour is al- 



most identical: only the higher temperature among 
and Tp plays a role, in consistence with the theoretical 
discussion. Therefore we compare the experimental re- 
sults at the base phonon temperature (of about 50 mK) 
to the theoretical results for Tp ^ Te in what follows. 

We studied Pep in the normal state as well by apply- 
ing a magnetic field of about 120 mT to suppress the 
superconductivity and measuring the partial Coulomb 
blockade (CB) signal Like in the superconducting 

state, two regimes are possible. In equilibrium the re- 
sults of Ref. apply. Under injection, the typical 
situation is such that Text ^ which we discuss now 
in more detail. The tunnelling rates in a state with an 
extra charge n for adding (-I-) or removing (— ) an elec- 
tron to the normal island with electrostatic energy change 
AF±(7i) = ±2Ec{n ± 1/2) T eV/2 are 

r^W = ;2]^ J dEh{E)[l-ME-AF^{n))]. (3) 

Here Ec = e^/2Cs is the charging energy of the island 
with the total capacitance Cs, and /i and /2 are the 
distributions on the source and target electrodes. For 
equilibrium distribution fi{E) = (1 -|- gS/fcaTi^-i ^[f^T^ 
Ti = T2, Eq. dSI) yields the result of Ref. Here we 

have the opposite limit of low bath temperature. Text = 
Ti < T2 = Te. For Ti = 0, fi{E) = l~e{E), yielding 
r±(?i) = (/cbT/c^Pt) ln(l + e-^-P*(")/'=^«^<=). The cur- 
rent into the island is / = e^^ ^^ (T(7i)[r+(?i) — r^(n)] 
where cr(rt) is the probability of having n extra electrons 
on the island. Since J2'^=-oo n-'^in-) = by symmetry, 
and X]^^L-oo ^(") ~ 1' fi'^'i ^'^^ the differential con- 
ductance up to the first order in Ec/kBTg 



Gt 2kBT,cosh\eV/4kBT,)' ^ ' 

The depth of the conductance minimum at 1/ = is 
AG/Gt = Ec/2kBT^ which is 50% larger than that 
in the equal temperature case. To find the width at 
half minimum we need to solve cosh^(ey±/4fcBTe) = 2. 
The fuU width is V^J^ = \V+-V-\ov V^J^ = 41n(3 + 

2\/2)/cBTe/e. This is about 65% of the equal-temperature 
value, V^f^ ~ 10.88/cBTe/e [H. 

Figure [3] is a collection of the data at the base temper- 
ature (~ 50 mK), in form of island temperature T^/Tc as 
a function of injected power. The superconducting state 
was measured for the three samples. The power has been 
normalized by that at Tc, to present data from different 
samples on the same footing. For samples A, B and C, 
P{Tc) = 14 nW, 3 nW, and 3 nW, respectively. The 
data on the three samples are mutually consistent. The 
quasiclassical result for a superconductor is shown by a 
solid line. The normal state data were taken for Sample 
C which is ideal for a measurement of the island temper- 
ature via partial CB: it has Ec/kB — 20 mK yielding an 
approximately 8% deep Coulomb dip in conductance at 
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INJECTED POWER TO THE ISLAND P/P(T^) 



FIG. 3: (color online) Energy relaxation from theory and ex- 
periment. The data in the superconducting state are from 
Sample A (squares), B (diamonds), and C (circles). The open 
triangles are from Sample C in the normal state. The solid 
line is the result of Eq. ([2} in the superconducting state. The 
dotted line indicates P/P{Tc) = (T/Tcf, and the dashed 
line P/P{Tc) = {T/Tcf. The inset shows three Coulomb 
peaks measured in the normal state under different levels of 
power injection: the solid lines are theoretical fits to them. 



zero bias at 90 mK (see the inset of Fig. [3]). The two 
large junctions were used for probing and the small ones 
for power injection. We first checked that the value V^^^ 
yields a good quantitative agreement with the equilib- 
rium temperature data over the whole range of the ex- 
periment. Next we measured the quasi-equilibrium elec- 
tronic temperature under injection. The low base tem- 
perature permits the use of the expression of V^J'^ above 
to extract in the range displayed in Fig. [3l Power law 
type behavior can be observed over the whole tempera- 
ture range Q.3Tc < Tg < Tq- The data approach those of 
the superconducting state near Tc — 1.45 K, as expected. 
The power law for Pop is, however, better approximated 
by (dashed line) instead of (dotted line) of Eq. 
([3]) , yielding a deviation of the same sign with respect to 
the basic theories as in the superconducting state. 

The experimental data demonstrate that e-p coupling 
in a superconductor is weaker than in the normal state, 
by two orders of magnitude at Tg/Tc = 0.3. But, like 
in the relaxation time experiments in a superconductor 
0i 0, Si , the energy flux is larger than that from the qua- 
siclassical theory [ll [l^ . This observation could suggest 
that the electron relaxation rate both in the supercon- 
ducting and in the normal state might be sensitive to 
the microscopic quality and the impurity content of the 
particular film J^4]. The impurity effects on the e-p re- 
laxation are controlled by the parameter q£ where i is the 



electronic mean free path and q = kBTe/hu is the wave 
vector of an emitted phonon with energy of the order 
of the electronic temperature. With the speed of sound 
u ^ 5000 m/s and £ 20 nm in our samples, we have 
q£ ~ 0.5 Tc. Thus the impurity effects can become 
essential below 1 K. More theoretical studies are thus 
needed of the impurity effects on the e-p interaction in 
the superconducting state. 

Our experiments were performed on three samples with 
very different lengths and junction parameters but they 
yielded essentially identical results when normalized by 
the island volume. Therefore we believe that issues like 
thermal gradients, slow electron-electron relaxation, and 
the presence of tunnel contacts have only a minor influ- 
ence on the results. The data thus yield the intrinsic en- 
ergy relaxation of QP:s in the superconducting and in the 
normal state. In summary, the experiment follows quali- 
tatively the theoretical model that we presented. On the 
quantitative level there is a substantial discrepancy espe- 
cially for superconductors, which would imply that one 
needs to invoke an extra relaxation channel to account 
for. Solution of this quantitative disagreement and ex- 
periments at still lower temperatures remain as topics of 
future work. 
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